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Abstract 

j Large assemblings in space require the ability to manage Rendezvous and Docking operations. In future these techniques 
1 will be required for the gradual build up of big telecommunication platforms in the geostationary orbit. 

i The paper discusses the use of Fuzzy Logic to model and implement a control system for the docking/berthing of two sat- 
| ellites in geostationary orbit. The system mounted in a chaser vehicle determines the actual state of both satellites and gen- 
t erates torques to execute maneuvers to establish the structural latching. The paper describes the proximity operations to 
1 collocate the two satellites in the same orbital window, the Fuzzy guidance and navigation of the chaser approaching the 
* target and the final Fuzzy berthing. 

I The Fuzzy Logic system represents a Knowledge Based Controller that realizes the close loop operations autonomously 
replacing the conventional control algorithms. The goal is to produce smooth control actions in the proximity of the target 
and during the docking to avoid disturbance torques in the final assembly orbit. 

The knowledge of the Fuzzy controller consists of a data base of rules and the definitions of the fuzzy sets. The knowledge 
of an experienced spacecraft controller is captured into a set of rules forming the Rules Data Base. 
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Introduction 

Since several years ago the number of spacecrafts in 
the geostationary orbit has increased considerably. 
Large satellites with longer life times are now the 
trends in the market. A big part of this development is 
due to commercial telecommunications companies 
with a growth on circuits demand of approximate 
10% per year. 

For the time being the satisfaction of this demand is 
fulfilled by a progressive tightening of the East-West 
deadband and the construction of complex clusters 


with collision avoidance strategies. But the concept of 
satellite cluster is only attractive if the functioning of 
the set of spacecrafts in the same orbital window is 
seen by the user as only one payload. 

Current practice for orbital windows in GEO assumes 
a square of about ±0.1° in latitude and longitude but 
this is gradually being shrunk to ±0.05°. This demand 
in orbital space has a big impact on spacecraft de- 
sign, ground station design and station keeping op- 
erations; the spacecraft must incorporate more effi- 
cient propulsion systems and satellite-to-satellite 
tracking devices (ESA 10035), the ground station has 
to increase in complexity to allow higher accuracy in 
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orbit determination and operations have to become 
more complex for inclination and eccentricity mainte- 
nance. Some solutions arise: satellite clusters (with 
heavy workload in ground operations or autonomous 
station keeping) or the gradual assembling of big plat- 
forms in orbit. The first technique requires orbit de- 
termination methods cm board the satellite (where 
ground operations workload is considerably reduced) 
or on ground. The second requires the ability to mas- 
ter rendezvous and docking operations but once the 
assembly is formed the control remains purely con- 
ventional. 

The advantages of a big platform formed by the join- 
ing of several pieces instead of a “dancing” cluster of 
satellites can be listed as follows: 

1 . reduction to zero the risk of collisions. 

2. reduction in ground station workload. 

3. possibility to add as many pieces as desired in- 
creasing by far the capacity of the initial pay- 
load. 

On the other hand the satellites must to be constructed 
with docking mechanisms and autonomous control 
systems. 

Rendezvous Mission Design 

Among the main three techniques used for a chaser 
approaching a target (R-bar, V-bar, inertial) the V- 
bar is the most common used (Leonard 89) for its 
stability properties. With this technique less propel- 
lant than R-bar is required but plume impingement 
problems appear. Here only V-bar option will be 
considered. For a mission of this type the following 
assumptions are made: 

• the targ et: it is a geostationary communications 
satellite. The target is passive and has the neces- 
sary mechanical-electrical elements for the dock- 
ing of the chaser. It is maintained in the centre of 
its orbital window within a specified inclination 
and eccentricity. The target is three-axis stabilized 
via a double-gimbaled, bias momentum control 
system during the complete rendezvous mission. 

• the chaser , it is another geostationary communi- 
cations satellite with identical shape and mass. 
The chaser has a control system that allows the 
rendezvous and soft docking with the target in an 
autonomous form (no man-in-the-loop capabilities 


are considered). If the rendezvous fails the chaser 
returns to a safety position. The chaser is ap- 
proaching the target using the V-bar technique 
where the docking axis is along the velocity vec- 
tor. 

• the environment : the chosen reference coordinate 
system is the Local Vertical, Local Horizontal 
(LVLH) (fig. 1). That is: +X in the direction of 
target flight, +Z in the direction of center of Earth 
and +Y towards the Earth south pole. The orbital 
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Figure 1. Final approach and injection into orbital window 

window measures are ±0.1° in longitude and 
latitude. This window is maintained by means of 
ground operations with a full antenna coverage in 
all occasions. The control centre could be located 
in ESOC (Germany) having an antenna which 
monitors the first part of the operation. 

Under these assumptions the type of mission design 
for an assembling in GEO of two satellites will not be 
different from the conventional design during the 
launch, transfer and drift orbits and the station ac- 
quisition phase (Pocha 87). The difference resides in 
the guidance and navigation of the chaser approach- 
ing die target and the final docking. 

Launch and GTO Orbit 

The launch can be performed from Kourou Space 
Centre in French Guyana using an Ariane 4 expend- 
able three-stage rocket. Using this rocket the inclina- 
tion of the transfer orbit is 7° . At the apogee of the 
transfer orbit the apogee motor is fired acquiring a 
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near geosynchronous orbit (the drift orbit). 

Drift Orbit and Station Acquisition 
During the drift orbit the chaser acquires its opera- 
tional three-axis stabilized attitude pointing involving 
de-spin. Sun acquisition and Earth acquisition. For 
this orbit a set of 4 ESA ground stations to track the 
target can be assumed. 

Once the satellite has acquired its operational attitude 
stabilization it must be placed closed to the target 
orbital longitude and with a specific orbit inclination. 
This operation is accomplished using ground tracking 
support. In this case only one station is involved. 

Rendezvous & Docking Operations 

These are defined as the set of operations to close up 

and dock two spacecrafts. In a typical mission they 



Figure 2. Rendezvous sequence timeline 


comprise (fig 2): homing, final approach, docking and 
structural latching (Wohlke 92, Pairot 92). 

Homing 

This phase starts with the target presence acquisition 
by the large range sensor (S-band radar) mounted in 
the chaser. That happens approximately at -100 Km 
(LVLH system) behind and slightly below the target. 
In this phase although the ground station supplies 
navigation data, the radar works in back-up mode 
giving R, A z and El of the target. During this phase 
only translation movement of both target and chaser 
is considered. 

At -60 Km the chaser crosses the orbital window and 


in that moment the medium range sensor (laser) 
locks -cxi the target. During this part the laser sensor 
and radar sensor are nominal. Now the attitude 
pointing of the chaser starts to be important. 

The homing phase extends to a distance of -1 Km 
(behind the target). In that moment another phase 
starts. 

Final Approach 

This phase comprises the close up of the chaser from 
-1 Km to -1 m. During this time the short range sen- 
sor (camera) localizes a specific mark in the target 
(Ho 93). In this phase the camera functions in nomi- 
nal mode whereas the laser and the radar are in 
backup mode. 

Docking 

This phase starts at -1 m from the target and ends just 
a few centimeters from it, before the latching. Four 
close up sensors mounted in a cruciform way in the 
chaser west platform side allow the fine docking. 

Latching 

Four latches mounted clockwise to the close up sen- 
sors will fit into four handles that close when the 
proximity operations are finished (Fehse 85). In this 
phase the sensors and the camera operate in parallel. 


Latch y 



Figure 3 . Docking mechanism 


The laser and radar data are not considered (fig. 3). 
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The physics of the system 

The translational motion of a spacecraft system in 
circular orbit can be described using the Clohessy- 
Wiltshire equations. Those are linear differential 
equations with time constant coefficients which de- 
scribe the movement of two masses in a circular orbit 
around a third object. They were programmed in the 
rendezvous guidance computer used in the Gemini 
mission (1962) and still provide short-range maneu- 
ver computation for the Shuttle. Nowadays the range 
of validity of these equations have bear extended by 
introducing special coelliptic coordinates to generalize 
the LVLH reference system. 


The Clohessy-Wiltshire equations have as input the 
initial position and velocity of the chaser. The output 
is the position and velocity of the chaser after a time 
interval. To apply these equations the two rotating 
bodies must have a small mass in comparison to the 
non-rotating body. In addition the target is located in 
the origin of the rotating coordinate system. The po- 
sition and velocity of the chaser are given in this 
LVLH reference system. 
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where r\j 2 are the distances of target and chaser from 
the Earth centre, }l is the product G*M«nh and f is 
the perturbing force. This equation linearised and 
written in terms of Cartesian coordinates (LVLH 
system) gives: 


x-2 COz = f x 


>+a> 2 }= f. 


2 + 2(0 x- 3 co 2 z = f z 

These equations are not solvable in general but in 
some special cases it is possible to derive an analyti- 
cal solution. If the external forces are zero (f = 0) the 
equations can be applied to a rendezvous in high LEO 
orbits (Shuttle, Soyuz-MIR) giving excellent results. 
(Malyshev 94, Brown 94). For a rendezvous in low 
LEO the force f is the atmospheric drag, etc. 


In the geostationary case several forces (apart from 
the Earth gravity) perturb the orbit of satellites: the 
Moon, Sun, Earth triaxiality, solar wind, etc. The 
Moon-Sun combined effect causes the orbit inclina- 


tion to grow 0.85° yr 1 in average. The Earth’s oblate- 
ness causes a precession of the ascending node of 
about 4.9° yz '. This effect is noticeable for non in- 
clined orbits but it is negligible for orbits with 0° in- 
clination. The Earth triaxiality causes a longitudinal 
acceleration towards GEO points at 79°E and 
107.6°W. An average value in acceleration of 0.001° 
per day 2 towards these points is typical. Finally the 
solar radiation pressure varies as the inverse square 
of the distance from the Sun. This force depends on 
die type of spacecraft surface and can be estimated 
as: 

Fgp =k* ;4 *1.6* cos/' 

The magnitude of this force along pitch axis is 10" 4 
•cos cOgeo *t Nm and the magnitude along roll and yaw 
axis is 10‘ 5 .cos ©goo.tNm. 


Once the chaser is close to the target during the final 
approach phase the force of the Earth gravity can be 
considered uniform. In this case the reference system 
becomes inertial body reference system and the equa- 
tions become: 
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where again f is the perturbing force and x, y, z are 



Figure 4, Thrust location in chaser 
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the chaser coordinates with respect to the LVLH cen- 
tre. The target system physics of movement can be 
described using body coordinate equations as well. 

The translation of both vehicles is realized by means 
of thrust impulses. The thruster system consists of 8 
pairs localized at four surfaces as shown in (fig. 4). 
Pairs (5,6) will impulse the satellite in the -X direc- 
tion and pairs (7,8) in +X direction, etc. 


Direction of Eorih and cho w Movomont 



Figure 5. Final closing body coordinate system 


The attitude control of both chaser and target is de- 
signed as a three axis stabilized, double gimbaled, 
bias momemtum control system (Kaplan 76, Wertz 
78). The corresponding Euler equation is 

(F ' +F »H£l +w * xh ‘ 

where S’ is a fixed reference system with origin in the 
centre of satellite mass and S is a reference system 
with the same origin that rotates with the satellite. S’ 
and S are defined as the LVLH system. F* is the dis- 
turbance force (solar pressure, thruster misalignment, 
etc.) and F„ is the force due to gravity gradient. w„ 
is the angular velocity of S respect to S' and h is the 
total angular momentum of the spacecraft (body + 
gimbaled wheel). 


Pitch equation is simple and decoupled from yaw and 
roll assuming a symmetric satellite 

but yaw and roll equations are coupled: 

F a =l x 4>+(O ss hJ+h w \i/+h x -CQ ss h t 

F*=l,V+(o ss h w \i/-h w <l>+h,-w ss h I 

where /» l y 1, are the principal moments of inertia 

and the control torques (h x ,h > ,h t ) are produced 
through gimbals defections (fig. 5). 

Fuzzy Control: System Model 

Figure 6 shows the proposed Guidance, Navigation 
and Control System (GNC). The navigation block 
calculates the actual state of both spacecrafts and 
their relative measurements. The guidance and control 
blocks are together in the figure. The guidance part 
calculates the future state of the spacecraft to achieve 



Figure 6. Guidance, navigation and control system 


the desired trajectory and the control part calculates 
the desired control torques to achieve this trajectory. 

It must be noticed that the guidance and control 
(G&C) system is divided in two parts: the G&C of 
the translation movement and the G&C of the rota- 
tional movement. For this study both movements are 
supposed to be de-coupled and the influence of one on 
the other is considered negligible. However we will 
see that the attitude motion influence on the G&C 
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part of the translation motion. 

Controller Implementation 

The fuzzy controller (FC) must substitute the guid- 
ance and control parts. The goal of the FC is to be 
able to generate the control signals derived from the 
sensor measurements (Lee 90). The steps involved in 
the construction of this FC are the following: 

1. define input and output variables. 

2. define universe of discourse for all variables. 

3. define the Rule Data Base (RDB). 

4. define the Inference Engine 

Once the FC is constructed it must be coded and in- 
cluded in the simulation. 

For all the four previous points a fundamental step in 
the implementation of the FC is to capture the knowl- 
edge of an expert. This knowledge must be applied to 
guide and to control the vehicle. The expert will help 
defining the way to determine the state of the target- 
chaser system and the way to use the commands for 
each situation. The inference engine can be an ap- 
proximate reasoning kernel based on already pro- 
posed systems (Buckley 92, King 94). 

Construction of IFS and OFS 

In this case the input variables are the sensor meas- 
urements (positions, velocities, etc.) and the output 
variables are the firing of thrust (thrust position and 
time of fire) and the attitude angles and rates. The 
input variables will be represented in the Input Fuzzy 
Sets (IFS) system and the output variables will be 
represented in the Output Fuzzy Set (OFS) system 
(Drianov 93). 

For the translation movement it is simpler to use input 
variables in polar coordinates rather than Cartesian: 
an aircraft pilot is able to measure azimut and eleva- 
tion angles and rates and distances and distances 
rates. The rotational movement can be represented by 
the attitude angles and angles rates. The output for a 
human controller could be the deflection of a joystick 
to fire a particular thruster and the defection time for 
acceleration controlled devices or the amount of 
pulses in pulse control devices. 

Defining the control variables 

For the translation movement of the chaser the input 
variables of the FC system will be azimut (a), eleva- 


tion (e), azimut rate (a), elevation rate (e), range 
(r) and range rate (r ). For the rotational part the in- 
put variables will be pitch (0), pitch rate ( 0 ), roll 

(ij>), roll rate (<j>), yaw (y) and yaw rate (y ). The 
output variables will be the amount of firing time (ft) 
and the position of fired thruster (pt). 

The universes of discourse of these variables are as 
follows: 

a e [-*,*] in rad, a e [-2,2] in °/sec 
e e [-*,*] in rad, e e [-2,2] in °/sec, 
r e [O.ILnJin Km, r e [-10,10] in m/s 
0 e [-*,*] in rad, 0 e f-0.5,0.5] in °/sec 
e [-w,x] in rad, <}> € [-0.5, 0.5] in °/sec 

ye [-w,x] in rad, \|f e [-0.5, 0.5] in °/sec 
Bum e [-10,10] in dcm/sec 

Defining the Fuzzy Sets 

The fuzzy subdivisions of the universe of discourse of 
any variable are called the fuzzy sets. Every ‘partial’ 
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Figure 7 . Azimut fuzzy set selection 


value of each variable can belong to more than one 
set. 


148 


For this study the following sets are assumed: 




* for azimut, azimut rate, elevation, elevation rate 
and range rate the fuzzy sets are Small Negative 
(SN), Small Positive (SP), Large Negative (LN) 
and Large Positive (LP) (fig. 7). 

* fix range the possible fuzzy sets are Small (S) and 
Large (L) as no more distinctions are needed. 

* for roll, pitch, yaw and its derivatives the fuzzy 
sets are Small Negative (SN), Small Positive 
(SP), Large Negative (LN) and Large Positive 
(LP) 



-h/2 -k/4 0 n/4 k/2 


Figure 8. Fuzzy representation for azimut 


Special attention must be put when creating the size 
and type for the fuzzy membership functions. De- 
pending on the position along the universe of dis- 
course axis, the control actions will be more or less 
smooth. A rule for relative good control states that the 
functions should overlap 25%. The profile of the 
function can be trapezoid, bell shaped, gaussian, etc. 
A first approach in the definition of the membership 
functions for the problem of the rendezvous is shown 
in figure 8. 


=> the control of each axis is carried out in an inde- 
pendent manner: firings in Z direction control azi- 
mut and firings in Y direction control elevation. 
The firings over each axis are calculated taking 
into account angles and its derivatives. 

=> the control in azimut and elevation is not symmet- 
ric: due to the Clohessy-Wiltshire equations of 
movement firings in X direction will ‘elevate’ the 
chaser in its path towards the target. It is neces- 
sary to realize compensation firings in the +Z di- 
rection. However, the control in elevation is sym- 
metric with respect to the +X axis. 

=> the tendency of human controllers is to make azi- 
mut and elevation equal to 0° during the final ap- 
proach phase (typically at -10 or -20 m of the tar- 
get) to be able to reduce control workload and fo- 
cus <Ni fine range rate control. This scheme com- 
plies with an intuitive proportional navigation 
guidance towards the point of starting axis trans- 
lation. 

=> the amount and size of the firings depend propor- 
tionally on the distance: different firing strategies 
must be followed depending on the distance to the 
target. When the chaser is in homing phase the 
amount of firings is low (typically 1 every 5 min- 
utes) and the size is big. In final approach phase 
the firings amount starts to increase (1 every mi- 
nute) but the size gets small. During docking the 
firings are frequent with very small size. 

=> attitude angles and angles rates are controlled 
during the translation along axis phase. 


The azimut of the chaser (movement along X axis) is Constructing the knowledge data base 

subdivided in the four fuzzy sets (LN, SN, SP, LP). The rules data base are given in the following form: 

For small angles (a e [-tc/4,ic/4]) the exponential if ANTECEDENT then CONSEQUENT 

gaussian shapes overlap a little bit more than 25 % where ANTECEDENT and CONSEQUENT are any corn- 

trying to smooth thruster actions within these limits. position of statements (and, or, not, etc). The state- 

ments contain declarations of associations of fuzzy 
Rules D a t a Base variables to fuzzy sets (fig. 9). 

The rules data base is the kernel of the knowledge 

base controller (KBC). Thanks to this data base the These rules are grouped in two categories and every 

FC will incorporate an experience which can wily be category contains several groups of control. One 

realized in the corresponding analytic model by means category corresponds to a long distance between tar- 

of manual operations. In this case the KBC imple- get and chaser and the other corresponds to short 

ments the close loop control actions substituting the distance between target and chaser. This distinction 

operator (Tolk 94).. allows to fire different sets of rules belonging to 

translation or rotational movement: during long dis- 
Capturing the knowledge tances (r is large) no rule of rotational movement is 

Here the expert tells us the following: fired. When the chaser is getting closer to the target 
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the relative attitude of both vehicles starts to be im- 



-2 0 2 

and 

azimut rate is LP 


Figure 9. Azimut rule 

portant. 

For each category the rules are grouped depending on 
the control action they generate: for the translational 
movement there are rules for controlling azimut, rules 
for controlling elevation and rules for controlling 
distance; for the rotational movement there are rules 
for controlling pitch, rules for controlling roll and 
rules for controlling yaw; 


Simulations 

For the simulations, the size of the satellites is as- 
sumed to be about 15 m high from the extreme ends 
of the solar panels and the body size is about 3 m 3 
(e g. ESA ECS satellite). The mass is 1000 Kgr. The 
matrix of moments of inertia is diagonal. The 
thruster misalignment torque is of an order of magni- 
tude less than the solar pressure. The attitude pointing 
requirements for the soft docking are 0.005° in pitch 
and roll and 0.4° in yaw. The docking velocity is 1 
cm/s. 

The plant description is simulated using the Clohessy- 
Wiltshire equations for the chaser influenced by a 
constant force in the -Y direction (Moon-Sun attrac- 
tion) during the first part of the rendezvous. During 
the final approach phase the reference system is iner- 
tial body and the equations get linear. The target 
translation is simulated with a constant force in the 
same direction. The attitude dynamics of both space- 


crafts are simulated using the Euler equations with 
gravity gradient and solar pressure perturbations. 

The FC controller is implemented in the following 
way: 

• The control variables together with their universes 
of discourse are as described before. 

• The shape of the membership functions is an ex- 
ponential gaussian function. 

• The rule data base is composed of 32 rules (20 for 
translation and 12 for rotation). 

• The inference engine is programmed using the 
Mamdani's Min-Max mechanism (Mamdani 74); 
the AND operator is chosen as the minimum of 
two weight antecedents instead of its multiplica- 
tion. 

• The defuzzyfication strategy used is the centre of 
gravity computation. 

The pseudocode of the simulation is as follows: 
initialize all; 

create rules for the rules data base; 
step = 1; 
bigloop 

for t = 1 to time to fire again 
draw positions & velocities; 
chaser kinematics; 
target kinematics; 
chaser attitude dynamics; 
target attitude dynamics; 
sensor measurements; 
store everything for drawing; 
end for; 

Fuzzy control computations; 
ifR = 0 then out 
else step = step+1; 
end big loop; 

Figure 10 shows a simulation run. The chaser was 
initially located in the coordinates (-1,0,1) behind the 
target at -lKm. The firing of thrusters try to close the 
chaser to the docking axis as soon as possible. This is 
practically achieved at -400 m. The translation along 
the axis occurred at -10 m. The azimut of the chaser 
decayed from 45° to 10° in a nearly exponential 
curve. From 10° to 0° the azimut waved up and 
down. This movement corresponds to the different 
trials to align with the docking axis. 

Figures 11a and lib show another simulation run. 
This time the chaser was initially located in the coor- 
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Figure 10. Simulation run (in-plane motion) 

dinates (-0.05,-0.05,0.05) behind the target at -500 m 
and below the target at -500 m as well. The firing of 
thrusters, for the in-plane motion (fig. 11a), try to 
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Figure 11a. Simulation run (azimut) in-plane motion 

close the chaser to the docking axis as soon as possi- 
ble. For the out-of-plane motion (fig. 1 lb) the firings 
were almost linear (because the angular velocity of 


the orbit is very low). The translation along the axis 
occurred at -1 1 m. The azimut of the chaser decayed 
again from 45° to 10° in a nearly exponential curve. 
From 10° to 0° the azimut waved up and down as 
well. The elevation of the chaser decayed from 45° to 
0° in a nearly linear curve. From 10° to 0° the eleva- 
tion did not waved up and down. 



x 

Figure Jib. Simulation run (elevation) out-of-plane motion 


Simulations were developed on a PC compatible 
486/66DX2 portable type computer. The code was 
written in MATLAB 4.0® for Microsoft Windows™. 
The fuzzy inference engine with approximate reason- 
ing was implemented using FISMAT, the Fuzzy In- 
ference Systems toolboox for MATLAB developed 
by A. Lofti. FISMAT provided fuzzy logic opera- 
tions and different methods of approximate reasoning. 

Contingencies 

Contingency situations can happen when the docking 
velocity is too high, the docking axis is not in perfect 
alignment with the V-bar axis or when the attitude of 
both satellites differ in an angle bigger than the one 
specified for the latching mechanism. 

Easily, the FC can be extended to handle these situa- 
tions. For an hazardous impact the thrust pair (1,3) 
must be used to produce an exponential breaking 
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during the translation along axis. This situation was 
taken into account with the rule 

if range rate is L then Bum (1,3) is SN 
For the situation with the docking axis not in perfect 
alignment with the V-bar axis a pseudo proportional 
navigation towards the starting point of the axis 
translation was implemented with the rule 

if a is LP and a rate is IP then Bum (5, 6) is IP 
Finally when attitude angles are out of range of the 
specifications of the latching mechanism some cor- 
rective measures can be taken during the translation 
along-axis phase. 

Conclusions 

Up to date, servicing is greatly dependent of human 
operations in space. Therefore, the on-orbit accessi- 
bility factors of the satellite to be serviced is a pri- 
mary consideration. Support capabilities for human 
presence in satellite maintenance tasks in GEO will 
not be widely available until the first decade of the 
21st century (Waltz 90). In the meantime autonomous 
rendezvous and docking arise as a good alternative. 

Fuzzy logic emulates the behavior of human opera- 
tors for complex control tasks. A fuzzy logic control- 
ler embedded in a guidance, navigation and control 
system of a spacecraft can realize autonomously the 
close loop operations replacing the conventional crisp 
control algorithms. The fuzzy controller produces soft 
control actions in the proximity of the target and 
during the docking to avoid disturbance torques in the 
final assembly. 

Fuzzy controllers can be programmed, tested and 
qualified for flight (Daley 85, 87). Its rule data base 
can be constructed with the help of an expert and re- 
fined in simulations. Fuzzy controllers are easily re- 
configurabled for different type of missions and their 
performance is robust to changes. 
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